arXiv:1504.04096vl [hep-lat] 16 Apr 2015 


PROCEEDINGS 

OF SCIENCE 



Canonical approach to the finite density QCD with 
winding number expansion 


Atsushi Nakamura 

RIISE, Hiroshima University 
Higashi-Hiroshima 739-8521, Japan 


E-mail: nakamura@riise . hiroshima-u .ac.jp 


Shotaro Oka 

Institute of Theoretical Physics, Department of Physics, Rikkyo University 
Toshima-ku, Tokyo 171-8501, Japan 
E-mail: |okasho-hato@rikkyo .ac.jp 


Yusuke Taniguchi* 

Graduate School of Pure and Applied Sciences, University ofTsukuba 
Tsukuba, Ibaraki 305-8571, Japan 


E-mail: taniqchi@het. ph ■ tsukuba . ac . jp 


The canonical partition function is related to the grand canonical one through the fugacity expan¬ 
sion and is known to have no sign problem. In this paper we perform the fugacity expansion by 
a method of the hopping parameter expansion in temporal direction for the lattice QCD: winding 
number expansion. The canonical partition function is constructed for Nf = 2 QCD starting from 
gauge configurations at zero chemical potential. After derivation of the canonical partition func¬ 
tion we calculate hadronic observables like chiral condensate and quark number density and the 
pressure at the real chemical potential. 


The 32nd International Symposium on Lattice Field Theory 
23-28 June, 2014 

Columbia University New York, NY 


* Speaker. 


(c) Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 


http://pos.sissa.it/ 













Canonical approach to the finite density QCD with winding number expansion 


Yusuke Taniguchi 


1. Introduction 

The grand canonical ensemble is a difficult subject to treat in lattice QCD because of the 
sign problem. The canonical partition function is related to the grand canonical one through the 
fugacity expansion and is known to have no sign problem as will be briefly reviewed in Sec. 2. 
In this paper we perform the fugacity expansion by a method of the hopping parameter expansion 
as is discussed in Sec. 3. The canonical partition function is evaluated for Nf = 2 QCD starting 
from gauge configurations at zero chemical potential in Sec. 5. After derivation of the canonical 
partition function we study the chemical potential dependences of hadronic observables like chiral 
condensate and quark number density in Sec. 6. 


2. Canonical partition function 

It is a well known fact that the grand canonical ensemble and the canonical one are equivalent 
each other. This is shown by a simple equation to relate the grand canonical partition function 
Zq(h,T,V) and the canonical Zc(n.T.V) 

oo 

Z G (fi,T,V) = £ Zc(n,T,V)S n , £ = e ^ T . (2.1) 

n——oo 

This is a so called fugacity expansion formula with the fugacity q. The inverse of this expansion is 
given by using the Cauchy’s integral theorem 

Z c (n : T,V) = j ^- n ~ l Z G (^T : V), (2.2) 

where we assume that the partition function Z ( ,(q. T. V) has singularities only at q = 0,oo corre¬ 
sponding to trivial points fi/T = ±°° and adopt an appropriate contour around the origin. Here we 
notice that a phase transition point is not a singularity of the partition function but it is rather a 
zero of Z G (%,T,V) in V —y °o limit (Lee-Yang zeros) and does not affect the Cauchy’s integral. 

Now it is free to change the contour to a unit circle q = e ,d and the contour integral turns out 
to be a Fourier transformation [[I]] 

r 27T fif) 

Z C (n,T,V) = / —e- md Z G {e' d J,V). (2.3) 

Jo 2% 

A standard Monte Carlo simulation is possible for the lattice QCD since the chemical potential is 
set to be pure imaginary fi/T =id and the Boltzmann weight is real positive for two flavors. 

However a difficulty of the sign problem is preserved unfortunately since there remains a 
frequent cancellation in plus and minus sign in a numerical Fourier transformation. In order to get 
the canonical partition function accurately for large quark number n we need a very fine sampling 
of Z G (e ,e ,T,V) in 0. This requires a heavy computational cost because we need to evaluate the 
Dirac operator determinant for the QCD grand partition function. In this paper we shall solve this 
problem by a direct expansion of the Dirac determinant in terms of the fugacity. 
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3. Winding number expansion 


We consider the lattice QCD grand partition function given in the path integral form 

Z g {$,T,V) = I DUDetD w (^U)e~ s ^ u \ (3.1) 

where we adopt the Wilson Dirac operator 

3 

D w = 1-kQ, Q = J £{Qt + Qfi)+e^Qi + e-^'Q^, (3.2) 

i=i 

{Qp) nm = I 1 “ %) U p( n )8m,n+th {Qg) nm = (l + Yu) (3-3) 

Both the chemical potential e~‘ la and the hopping parameter K appear in front of the temporal 
hopping term simultaneously. The fugacity expansion of the determinant shall be given by using 
the hopping parameter expansion. 

Instead of the determinant we consider the hopping parameter expansion of 


TrLog£%(jU) = TrLog (/ - kQ) 



(3.4) 


Since every term has a trace the quark hoppings should make a loop. A non-trivial chemical 
potential dependence appears when the quark hoppings make a loop in the temporal direction. If 
one of the term has a n times winding loop in positive temporal direction the chemical potential 
dependence is e^ anN ’ where N, is a temporal lattice length and this is nothing but e'— y = q". 
Counting the winding number in temporal direction for each term in ( |3~| ) we get the winding 
number expansion of the TrLog 


oo 

TrLog D w (ii)= £ w n q n 

Yl ——°° 


(3.5) 


Regrouping the summation we have a fugacity expansion of the determinant [|2|] 


DetDw(^;I7) = exp 



oo 

= E 

Yl ——oo 


(3-6) 


In this approach the fugacity dependence of the determinant or the partition function is given ana¬ 


lytically. A numerical Fourier transformation (|2.3|) can be executed safely at any high precision we 


want. In this paper the regrouping (p.6| ) is done by the Fourier transformation. 

The winding number expansion is done for gauge configurations generated at = 0 or purely 
imaginary value. This procedure corresponds to the standard reweighting method 




(3.7) 
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Table 1: (/3, k) for the numerical simulation. The PCAC Figure 1: The Polyakov loop as a function 
mass and m n /m p are also given. of /). 


4. Numerical setup 

We adopt the Iwasaki gauge action and the improved Wilson fermion action with the clover 
term. The number of flavors is set to Nf = 2 with degenerate masses. The APE stout smeared gauge 
link is used for those in the fermion action including the clover term. The number of smearing is 
four and the parameter is set to p = 0.1. The clover coefficient is fixed to csw = 11 [01- We adopt 
8 3 x 4 lattice. The numerical parameters (/3. k) are given in table 0. A wide range of /3 is covered 
from high temperature /3 =2.1 to low temperature 0.9. It seems to be that both the confining 
and deconfining region are well covered, which is inferred from a behavior of the Polyakov loop 
in Fig. | The hopping parameter is selected in order that the hopping parameter expansion works 
well, for which m K /m p turns out to be 0.7 - 0.9. Maximal order of the hopping parameter expansion 
is set to 480 so that max winding number in temporal direction is 120. 


5. Canonical and grand canonical partition function 


The first numerical result we get is the canonical partition function Zc{n,T,V). We plot 
log | Zc(tiB, T,V)/Zc{ 0, T,V)\/(VT 3 ) and its phase as a function of the baryon number ///; in Fig. 0. 
The partition function decays very rapidly with ns and its behavior changes drastically between 
/3 = 1.7 (magenta) and 1.3 (dark-green), which may correspond to a phase transition. The phase 
of the partition function is consistent with zero within the error bar for /3 = 1.9 as is shown in the 
right panel. 

The plot range is fixed by using the d’Alembert’s convergence condition 


lim 

n B ^oo 


ZcjnB +1) g 
Zc(hb) 


< 1 , 


(5.1) 


1 which gives the convergence radius for the fugacity q. The data at JS = 1.1 and 2.1 are plotted 
in right panel of Fig. |0 for example. We cut our data at « max where a monotonic decrease stops 
indicated by vertical blue and orange lines in the figure. The horizontal lines show maximal values 
of the fugacity expected to be within the convergence radius. By taking log the line gives our 
applicable limit for the baryon chemical potential fis/T. For example we can discuss physics 
safely at —10 < fis/T < 10 for /3 = 1.1 and —4 < Hb/T < 4 for /3 = 2.1 with our method. 

'For negative baryon number we adopt |Zc(«b — 1) /Zcfiig )|. 
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Figure 2: The canonical partition function as a function of the baryon number n B . Left panel is the absolute 
value \og\Z c (n B )/Z c {0)\/(yT 3 ). From the top fi = 2.1 (orange), 1.9 (red), 1.7 (magenta), 1.5 (green), 1.3 
(dark green), 1.1 (blue), 0.9 (black). Right panel is a phase of Zc(n B ) at /3 = 1.9. 



Figure 3: Log plot of | Zc(n B + l)/Zc(n B )\ as a function of the baryon number. The color and p correspon¬ 
dence is the same as in Fig. 2 for left panel. Right panel is that at jS = 1.1 and 2.1 for example. 


The second physical quantity is the grand partition function. By taking summation for —n max < 
n B < « max in ( p~T| ) we get the grand partition function for the real chemical potential. We plot 
(log |Zg(jU, T, V)/Zg(0, T,V)\)/(VT 3 ) in Fig. [i| as a function of the quark chemical potential. Ac¬ 
cording to the statistical physics the logarithm of the grand partition function divided by the spatial 
volume is the pressure 

L = L\ og z G . (5.2) 

The quantity plotted in Fig. [I] is the pressure (P(fJL/T) — P(0))/T 4 in the grand canonical ensemble 
normalized at /./ = 0. 


6. Hadronic observables 


In our procedure the fugacity expansion is based on the hopping parameter expansion. It may 
be possible to expand any hadronic operators in terms of the fugacity. We consider a fugacity 
expansion of a numerator of some operator VEV 


(0) G (^T,V) 


Og^JX) 

Zg&JXY 


( 6 . 1 ) 
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Figure 4: The grand partition function (\og\Zc{p,T,V)/Zc(0,T,V)\)/V as a function of the quark chemical 
potential p/T. The color and /3 correspondence is the same as in Fig. 2. 
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Figure 5: Left panel is the chiral condensate (yy) /( VT 3 )) in the canonical ensemble as a function of the 
baryon number. Right panel is that in the grand canonical ensemble as a function of the quark chemical 
potential fl/T. The color and fi correspondence is the same as in Fig. 2. 


O g (Z,T,V)= DU(0) quark ^)DetD w ^-U)e~ SG{u) = £ O c (n,T,V)Z n , (6.2) 

n =—oo 


where (0) quark is a Wick contraction of a hadronic operator O in quark fields. For example we 
consider the chiral condensate. It is easy to expand its Wick contraction in terms of the hopping 
parameter 

- -«(£) - - ijha) - %r r(t = (6 31 


Counting the winding number in temporal direction we get the last equality. Multiplying the deter¬ 
minant contribution ( |3.6| ) we apply the same regrouping procedure in Sec. || and get the fugacity 
expansion (5.2). 

Once we have two coefficients Zc(n,T,V) and Oc(n. T. V) the VEVs of the operator with the 
canonical and the grand canonical ensemble are available. The canonical ensemble VEV is given 
by taking the ratio of two coefficient 


(0) c (n,T,V) 


Qdnjy) 

Zc{n,T,Vy 


(6.4) 


The results for the chiral condensate — (yy) c /(VT :, ) are given in the left panel of Fig. | as a 
function of the baryon number. A VEV in the grand canonical ensemble is given by taking fugacity 
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Figure 6: Quark number density — (i/F t/r) G as a function of the quark chemical potential. 


summation with real chemical potential 


(0) g ^,T,V) 


Z„=-n m3X Oc(n,T,V )g" 
L”-„ max z c (n,r,y)^ 


(6.5) 


The chiral condensate — (i//i//) G /(Fr 3 ) is given in the right panel of Fig. || as a function of the 
quark chemical potential /i /T. The condensate in the figure is a bare quantity without renormal¬ 
ization. Since we adopted the Wilson fermion we have an additive correction for (i//i//j, which is 
not subtracted in this paper. 

From the right panel of Fig. |]the chiral restoration phase transition at finite chemical potential 
seems to be seen. A relatively large value around fi/T = 0 decreases at large chemical potential 
and a fall occurs rapidly at a narrow region. The would-be transition parameter /T becomes 
larger for lower temperature. 

In Fig. |6| the quark number density (i// ; y/j f , /(FT 3 )) in the grand canonical ensemble is plotted 
as a function of the quark chemical potential \i/T. 


7. Conclusion 

In this paper we performed the fugacity expansion of the grand partition function by using 
the hopping parameter expansion. This procedure seems to be valid for baryon numbers around 
ns ~ 30 for 8 3 x 4 lattice. The method is also applied to the numerator of VEVs of hadronic 
operators. Taking summation we get a VEV at the real chemical potential. As an example we 
evaluate the chiral condensate and show a phase transition like behavior at high chemical potential 
for low temperature region. 

This work is done for Zn Collaboration. This work is supported in part by Grants-in-Aid of 
the Ministry of Education (Nos. 26610072. 24340054, 22540265). This work is in part based on 
Bridge++ code (http://suchix.kek.jp/bridge/Lattice-code/). 
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